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Calculating  ground  freezing  is  of  considerable  theoretical  and 
practical  interest;  thus,  knowledge  of  the  depth  of  freezing  is 
necessary  in  solving  numerous  problems  in  construction  and  agriculture, 
for  instance:  determining  the  depth  at  which  water-supply  and  sewage 
pipes  are  laid,  determining  depths  for  laying  the  foundations  of 
buildings,  the  initial  depths  for  digging  drains  when  carrying  out 
reclamation  projects,  etc.  Therefore,  during  the  last  20  years  a 
number  of  studies  have  appeared  devoted  to  this  problem. 

This  article  gives  an  analysis  of  some  of  the  existing  methods 
of  calculation;  in  this  process  we  have  selected  only  those  which 
at  the  present  time  should  be  recommended  for  engineering  and  whose 
area  of  application  of  each  method  is  indicated.  In  addition,  new 
computation  functions  have  been  derived  which  should  be  used  in  a 
number  of  cases  associated  with  a considerable  thermal  influx  from 
thawed  ground. 

I.  Some  Observations  on  Existing  Computation  Methods 

We  will  examine  theoretical  solutions  to  the  problem  of  ground 
freezing  as  given  in  the  works  of  G.  Gerber,  N.  N.  Petrunichev  and 
G.  S.  Shadrin,  I.  A.  Chamyy  and  B.  V.  Proskuryakov. 

V 

Greber  [1]  formulates  the  problem  which  was  previously  solved 
Stefan,  as  follows:1  "at  the  initial  moment  in  time  all  points  of  a 
humid  soil  have  temperature  T.  Beginning  at  a certain  moment  in  time, 
on  the  soil  surface  constant  temperature  6 is  always  maintained.  Freezing 
of  the  soil  takes  place  at  temperature  tQ.  In  order  to  freeze  one 

volumetric  unit  of  soil,  a quantity  of  heat  equal  to  w is  required. 

We  will  find  an  equation  for  the  temperature  field  and  the  speed 
with  which  the  freezing  boundary  moves"  (see  Figure  1). 

The  solution  has  the  following  form: 

(1) 

where  x is  the  depth  of  freezing, 
t is  time, 
q is  a coefficient. 

In  order  to  find  the  coefficient  q which  determines  the  freezing 
speed,  it  is  necessary  to  solve  the  following  equation: 


^e  designations  are  shown  in  accordance  with  those  adopted  in  this 
article. 
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where 


_ ,J  n» 

a-k  • 4,1  ,.  « *“  2 

(2) 

Vo-6)- 

(3) 

(4) 

In  equations  (1),  (2),  (3)  and  (4)  the  letters  X,  c,  y and  a designate 
the  thermal  conductivity,  heat  capacity,  specific  weight  and  temperature 
conductivity  of  the  ground;  index  "l"  refers  to  frozen  ground,  and 
index  "2"  refers  to  thawed  ground. 


The  following  should  be  noted 
about  the  above-cited  solution  to 
the  problem  of  ground  freezing: 

1.  The  solution  is  given  for 
conditions  which  considerably 
restrict  its  practical  application; 
in  particular,  the  solution  is  not 
applicable  in  the  presence  of  one 
of  the  following  conditions: 

a)  in  the  presence  of  a 
snow  covering, 

b)  if  during  the  freezing 
process  the  temperature  of  the  soil 
surface  changes, 

c)  in  the  presence  of  a 
filtration  flow. 

2.  For  practical  purposes  it  is  necessary  to  introduce  the  assumption 
that  the  temperature  of  the  soil  surface  is  equal  to  the  air  temperature. 

3.  The  performance  of  the  calculations  is  associated  with  cumbersome 
mathematical  computations  which  are  required  to  determine  coefficient  q. 

4.  Within  the  boundaries  of  the  stated  conditions  the  problem 
solution  is  analytical.  Therefore  this  method  should  be  used  wherever 
there  is  no  snow  covering  or  filtration  flow  and  where  the  change  in  air 
temperature  can  be  averaged  in  time  and  the  initial  temperature  of  the 
ground  can  be  averaged  in  depth,  and  these  factors  can  be  assumed  to 

be  constant  in  the  calculation. 


Figure  i. 


2In  his  derivation  Greber  made  an  error  which  leads  to  a plus  sign 
in  front  of  the  last  term  in  the  equation.  !'.’e  draw  the  reader's 
attention  to  this  fact  especially  because  this  error  has  been 
consistently  repeated  up  to  the  present  time  by  other  authors. 
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In  order  to  overcome  the  difficulty  associated  with  determining 
coefficient  q,  we  propose  to  make  use  of  nomograms  which  we  have 
devised  (see  Figures  2,  3 and  4). 

The  unknown  proportionality  coefficient  q is  presented  in  the  form 
of  two  terms: 


9-9i-9t» 

From  equation  (2)  it  follows  that: 


e-Ml 

a\Mt) 


(km‘k*  1 — ' 


(5) 

(6) 
(7) 


where 


From  these  two  equations,  we  find 


(8) 

(9) 

(10) 


9-  KdfljAf,.  (H) 

The  nomogram  in  Figure  2 relates  to  the  case  of  ground  freezing 
in  the  absence  of  heat  flow  from  the  thawed  ground,  i.e.,  when  T * tQ 

and  consequently  ^ * °*  The  nomogram  in  Figure  3 relates  to  the  case 

where  the  ground  thaws  under  the  action  of  the  heat  influx  from  the 
thawed  ground  in  the  absence  of  a thermal  flow  in  the  frozen  ground, 
i.e.,  when  0 » t^  and  q^  ■ 0.  The  nomogram  in  Figure  4 relates  to 

the  general  case,  i.e.,  when  there  are  thermal  flows  in  both  the 
frozen  and  thawed  ground,  when  0 i tQ  i T.  The  method  of  utilizing 

the  nomograms  in  Figures  2 and  3 and  the  principles  of  their 
construction  do  not  require  explanation.  They  should  be  used  when  q1 

is  considerably  greater  (or  less)  than  q^.  With  regard  to  the  nomogram 

in  Figure  4,  some  remarks  are  in  order. 

Figure  4-a  (the  attachment)  shows  two  families  of  curves:  one 
curve  is  plotted  according  to  equation  (6),  and  the  second  is  plotted 
according  to  equation  (7).  In  Figure  4-b  (attachment  at  the  end  of 
the  book)  there  are  also  two  families  of  curves:  one  is  plotted 
according  to  equation  (10),  and  the  second  according  to  equation  (11). 
The  value  of  the  unknown  coefficient  q is  determined  in  the  following 
fashion. 
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Figure  3-  Thawing  — Thermal  Influx  Only  From  Thawed  Ground: 

_ <u. 

t *** 

******  i cf  <'t  **  k2  <0,00l»  then  ^2  " k2’ 

' >.lW 

We  use  the  initial  data  to  find  the  values  k^,  k^,  and  a^/a^ 

We  will  set  an  arbitrary  value  of  Mj,  and  then  in  Figure  4-b  we  will 

determine  the  values  q and  M2;  then,  in  Figure  4-a  we  find  the  values: 

q^,  corresponding  to  and  corresponding  to  M2,  (Translator's  Note: 

one  word  obscured]  difference  should  be  equal  to  q;  when  these  values 
do  not  coincide,  it  is  necessary  to  select  another  value  for  Mj. 

e 

N.  N.  Petrunichev  and  G.  S.  Shadrin  solved  the  freezing  problem 
by  utilizing  the  method  of  finite  differences  [2].  Let  us  note  the 
following  about  this  solution: 

1.  this  method  can  be  applied  not  only  for  the  conditions  of  the 
problem  formulated  by  Greber,  but  also 

a)  in  the  presence  of  a snow-covering, 

b)  if  during  freezing  the  air  temperature  changes, 

c)  in  the  presence  of  a filtration  flow, 
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d)  with  any  initial  ground  temperature  distribution  with  regard 

to  depth. 

Even  the  presence  of  all  of  these  factors  together  does  very  little 
to  complicate  the  problem  solving  technique. 

2.  This  method  makes  it  possible  to  allow  for  the  difference  between 
the  air  temperature  and  the  ground  (or  snow)  surface  temperature. 

3.  This  method  makes  it  necessary  to  determine  a large  number 
of  intermediate  freezing  depth  values,  and  the  performance  of  the 
calculation  frequently  becomes  extremely  cumbersome.3 

4.  This  method  should  be  used  when,  due  to  partial  and  considerable 
changes  in  air  temperature  the  course  of  the  temperature  can  be  reflected 
only  by  a multi-step  (approximately  more  than  five  steps)  line  or  when 
the  initial  temperature  distribution  throughout  the  depth  of  the  ground 
becomes  non-uniform,  and  cannot  be  averaged. 

I.  A.  Charnyy  [3]  and  B.  V.  Proskuryakov  [4]  have  developed 
simplified  analytical  methods  for  solving  the  freezing  problem;  in  this 
process  the  two  authors  have  made  the  following  uniform  assumptions: 

1.  averaging  of  the  temperature  gradient  in  the  frozen  ground; 

2.  when  calculating  the  advance  of  the  freezing  boundary  the 
depth  of  freezing  x^  is  determined  which  occurs  in  the  absence  of  a 

thermal  influx  from  the  thawed  ground;  then  layer  x^,  which  thaws  under 

the  action  of  the  thermal  influx  from  the  thawed  ground,  is  subtracted, 
i.e. , 


* *i  - (12) 

The  basic  difference  in  performing  the  calculation  according  to  the 
method  of  I.  A.  Charnyy  and  D.  V.  Proskuryakov  arises  in  connection 
with  determining  x^.  The  heat  balance  equation  which  is  related  to  a 

vertical  section  of  ground  with  a cross-section  area  equal  to  unity, 
according  to  Charnyy  has  the  following  form: 

9f-e#-<?g-Qh,  (I*) 

and  according  to  Proskuryakov  it  has  the  form 


Q«-Q 


h* 


(14) 


3This  drawback  is  alleviated  by  utilizing  Petrunichev's  proposed  method 
of  performing  tabular  integration  of  the  Fourier  equation  instead  of  the 
graphic  method  which  is  used;  this  makes  the  performance  of  the  computation 
somewhat  faster  and  more  convenient. 
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where  Qf  is  the  quantity  of  heat  liberated  when  the  ground  within  the 
limit  of  the  section  freezes, 

Qq  is  the  quantity  of  heat  which  escapes  to  the  atmosphere, 

Qg  is  the  change  in  the  heat  content  of  the  frozen  ground, 
is  the  thermal  influx  from  the  thawed  ground. 

In  solving  the  equation,  Chamyy  obtained: 


x,  - 


y v —Yici  ~2  -- 


(15) 


Proskuryakov  found: 


x,~Y -f. 


where 


/*-/?-£ — A , 


where  h is  the  thickness  of  the  snow  covering. 


(16) 


(17) 


As  is  the  thermal  conductivity  of  the  snow, 

a is  the  coefficient  of  snow-air  thermal  emission, 

x is  the  initial  depth  of  freezing. 

0 

It  must  be  emphasized  that  the  expressions  for  Qq  are  not  identical 

in  equations  (13)  and  (14).  This  is  due  to  the  fact  that  according  to 
Chamyy  and  Greber,  0 is  the  temperature  of  the  ground  surface4;  according 
to  Proskuryakov  9 is  the  temperature  of  the  air  and  the  thermal  resistance 
of  the  snow  and  the  inequality  between  the  temperature  of  the  snow  surface 
and  that  of  the  air  are  taken  into  account. 

In  the  absence  of  a filtration  flow  Chamyy  and  Proskuryakov 
determined  [Translator's  Note:  word  obscured]  by  the  same  means,  using 
the  solution  from  the  theory  of  heat  transmission  to  the  problem  of  a 
loss  of  heat  by  a rod  which  is  limited  on  one  side  and  which  at  the 
initial  moment  of  time  has  a uniform  temperature  at  all  its  points  and 
when  there  is  no  heat  emission  from  its  lateral  surface;  on  this 
assumption  the  authors  obtain: 


4In  his  work  Chamyy  made  no  special  examination  of  ground  freezing, 
but  rather  "the  advance  of  a flat  division  boundary  between  two  phases"; 
therefore  it  is  natural  that  the  problem  of  what  influence  the  snow 
covering  had  was  not  stated  there. 
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(18) 


x.,»ki{V*  — V t0). 


For  the  case  of  a filtration  flow  Proskuryakov  found: 

xs  — ~ k.  (r  -To).  (19) 

where  y is  the  length  of  the  filtration  flow  path  under  the  frozen 
ground  to  the  point  where  the  depth  of  freezing  is  determined; 

v is  the  speed  of  the  filtration  flow's  motion. 

With  regard  to  these  two  assumptions,  we  should  point  out: 

1.  First  assumption,  averaging  of  the  temperature  gradient  in  the 

frozen  ground,  causes  an  increased  temperature  gradient  at  the  freezing 

boundary  and  a reduction  of  this  gradient  at  the  ground  surface;  the 

value  Q.  which  is  calculated  according  to  heat  balance  equation  (13) 
t » 

is  reduced,  but  according  to  equation  (14)  this  value  is  increased. 

As  a result  computation  formula  (15)  gives  understated  values  of 
x,  and  formula  (16)  gives  exaggerated  results;  it  is  known  only  that  the 
true  values  are  found  within  these  limits  which  also  limit  the  possible 
error.  As  calculations  show,  this  error  does  not  exceed  30%  for 
practical  purposes.  Taking  into  account  the  usual  imprecision  of 
initial  computation  data,  we  must  acknowledge  that  the  introduction 
of  this  assumption  is  completely  justified. 

2.  With  regard  to  the  second  assumption  which  is  expressed  by 
equation  (12),  it  is  obvious  that  it  causes  some  increase  in  the  influence 
of  the  thermal  flow  from  the  thawed  ground;  a qualitative  evaluation  of 
this  factor  is  difficult,  and  therefore  we  will  return  to  this  problem 
later  after  we  obtain  new  computation  equations  which  are  derived 
without  the  aid  of  this  assumption. 

With  regard  to  determining  the  freezing  of  ground  by  means  of  simplified 
analytical  methods,  the  following  should  be  said: 

1.  The  methods  can  be  applied  if  the  initial  ground  temperature 
distribution  can  be  averaged  in  depth  and  if  the  course  of  the  air 
temperature  is  replaced  by  a staggered  line  within  the  limits  of 
three-five  steps. 

2.  Formulas  (16)  and  (18)  make  it  possible  to  perform  calculations 
in  the  presence  of  a snow  layer  and  with  allowance  for  the  presence 

of  a difference  between  the  temperature  of  the  ground  surface  and  the 
temperature  of  the  air.  In  case  there  is  a filtration  flow,  formula 
(18)  should  be  replaced  by  formula  (19). 

3.  In  case  there  is  no  snow  covering  and  the  temperature  of  the  soil 
surface  is  equal  to  the  air  temperature,  it  is  possible  to  replace 
formula  (16)  with  formula  (15).  As  has  already  been  indicated,  computation 
according  to  formula  (16)  gives  somewhat  exaggerated  results,  and 
computation  according  to  formula  (15)  provides  correspondingly  reduced 
values  of  the  freezing  depth. 
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In  solving  the  problem  of  freezing  in  the  presence  of  a snow  covering 
and  with  allowance  for  a difference  between  the  air  temperature  and  the 
temperature  on  the  ground  surface,  as  is  done  by  Proskuryakov,  but  in 
contrast  to  him,  by  utilizing  equation  (13)  we  obtain: 


where 


T—T„ 


i A (x.\ — xj)  + <V  (Xj  — x0)  + R In 


(20) 


A- 

A' 


R 


>v  _ 
V.x  (0  - /,)  ’ 


JL 

2a. 


-h  2.4/9; 


(21) 

(22) 

(23) 


We  neglected  to  derive  formula  (20)  since  in  principle  it  is  similar 
to  the  course  of  the  derivations  carried  out  by  Proskuryakov  and 
Chamyy  when  they  obtain  formulas  (15)  and  (16);  let  us  simply  point 
out  that  for  us 


c.y.  (0 — t,) 

r /*  x,  .* 

2 

l(x+0y  * 1 .t-r/S 

dX  ~2(X-P)-' 

dx. 


(24) 


Like  equation  (15)  equation  (20)  provides  smaller  values  of  x^, 

but  in  contrast  to  equation  (15)  the  problem  is  solved  with  allowance 
for  the  difference  between  the  air  temperature  and  the  temperature 
on  the  ground  surface,  as  well  as  in  the  presence  of  a snow  covering. 

II.  Deriving  Refined  Computation  Formulas 

The  simplified  analytical  method  cannot  be  used  in  some  cases 
since  sometimes  it  is  not  adequately  precise.  If  we  proceed  on  the 
basis  of  heat  balance  equation  (14)  and  do  not  introduce  assumption 
(12),  then  the  formulas  obtained  should  precisely  allow  for  the  action 
of  the  heat  flow  from  the  thawed  ground.  Otherwise  all  of  the  advantages 
and  disadvantages  of  the  above-described  simplified  analytical  method 
will  also  be  noted  here. 

Let  us  examine  two  possible  cases.  First  case.  The  problem  can  be 
formulated  as  follows. 

At  the  initial  moment  of  time  the  moist  ground  has  temperature  T 
at  all  points.  The  surface  of  the  ground  is  covered  with  a layer  of 
snow  of  thickness  h.  Beginning  at  a certain  moment  in  time,  the  air 
temperature  is  kept  equal  to  0.  The  ground  freezes  at  temperature  tQ. 

A quantity  of  heat  equal  to  w is  required  to  freeze  a volumetric  unit 
of  ground.  The  physical  constants  of  the  ground  and  the  snow  and  the 
snow-air  thermal  emission  factor  are  given.  It  is  necessary  to  determine 
the  speed  with  which  the  freezing  boundary  moves.  Under  these  conditions 


the  intensity  of  the  thermal  influx  from  the  thawed  ground  q is 
inversely  proportional  to  time.  T 

The  heat  influx  term  from  the  frozen  ground  will  be: 

Qf  (25) 

where 


<1 


t 


{T  ~ >.) 

Vx  " ’ 


The  other  terms  of  heat  balance  have  the  following  form: 


(26) 


i 


Qo 


Qf  = wdx, 

“**  *+7*- 


(27) 


(28) 


__  f.KThe  latter  equation  allows  only  for  convection  heat  exchange 
" S"?T  surfaca;  however,  it  is  not  particularly  difficult  flso 
of  thkC  allowan*e  foj;  the  tenns  of  solar  radiation  R and  the  radiation 
V surface.S*  For  thi®  Purpose  it  is  necessary  in  all 

0 bvTf^t?rPUtatl°nSJ°ureplaCe  the  true  £«»Perature  of  the  air 
9 by  a fictitious  one  which  is  calculated  according  to  the  equation: 


%-«  + 


R-S 


eouat?onSnI?itULing  exPreasions  (25).  (26),  (27)  and  (28)  into 
equation  (14),  after  transformation  we  obtain: 


* + 

for  k2,  see  expression  (4) . 


.i  *i  ^ 

2 Vx  + 2 (x  + ?) 


(29) 


i,- 

r 

(30) 

We  will  introduce  a new  variable 

(31) 

then 

• T-r»(x+fy, 

(32) 

dr  «2r(x+/0*<//-  + 2r?(* +/*)dx. 

(33) 

-10- 


By  substituting  (32)  and  (33)  into  (29),  after  simple  transformations 
we  will  have 


dx 


*.-Hy 

i+V+V* 


dr. 


By  integrating  this  equation  with  allowance  for  the  fact  that 
4k3  - k*  < 0: 

/Vyi-kt-Vkl— 4kt  \ V i*. 

ln(x+fl-ln  ; 

VI  +*,r+*y* 


from  this  we  find: 


1 ,r  + P\m 

V(r^p){r^n)  'r+n  } 

(34) 

where 

(35) 

'-&(>*  V 

(36) 

m— , 

2^.-3- 

(37) 

The  integration  constant  is  determined  from  the  following  initial 
condition;  when  t ■ tQ,  we  have 

x-X,, 
r _ . 

from  this  we  find 


C - (x,  + 0)  / (r, + p)  (r0 + n)  • (J^) 


(38) 


By  substituting  the  expression  for  the  integration  constant  into 
equation  ([Translator's  Note:  number  obscured])  we  finally  find: 


- - -t  v ( '•+* -V"(  r-±>  V"  -ft  (39) 

TPl  y(r7p)(7|liy  \r.+p)  \r+n  ) p 


Second  case.  The  problem  conditions  are  the  same  as  in  the  first 
case,  but  the  intensity  of  the  heat  influx  from  the  thawed  ground  has 
a constant  value  like,  for  instance  [word  obscured]  by  means  of  the 
filtration  flow  from  the  lower  surface  of  the  frozen  ground,  as 


-11- 


7T-  (T_'>)- 


Then  heat  balance  equation  (14)  will  have  the  form 

dr  + wdx  + qt<lt^>  0, 

* — 


where 


A-W-tJ+q^. 


We  integrate  equation  (41): 


r— JLI_,..li!ij.i.|B(^TV)  + c.  (4: 

We  determine  the  integration  constant  from  the  initial  condition: 
when  t * Tq,  we  have  x » x^;  consequently 

C - r#^-  -BL  xtt -r  iv  7— ■ 1“  (-4  + <7t  vo>- 

?t  ?t 


Finally  we  find: 


(fill  — A)  a- ; <•  w 


By  utilizing  equations  (39)  and  (44)  and  equations  (16),  (18)  and 
(19)  of  the  simplified  analytical  method,  we  perform  numerous  calculations 
which  showed  that  the  introduction  of  computations  carried  out 
according  to  equations  (39)  and  (44)  refined  the  solution  by  no  more 
than 

’ " T,^T  ,C0K*  (45) 

In  actuality  the  refinement  is  less  significant,  but  the  latter 
can  only  be  shown  by  pairing  the  results  of  calculations  carried  out 
according  to  both  methods.  Consequently,  if  for  example  n > 20\,  then 
the  calculations  have  to  be  carried  out  according  to  equations  (39) 
or  (44),  although  the  refinement  may  also  prove  to  be  small. 

III.  Areas  of  Application  for  the  Computation  Methods 

The  selection  of  the  computation  method  which  should  be  used 
depends  on  the  specific  given  conditions  under  which  freezing  takes 
place  and  on  the  required  precision  of  the  problem  solution. 

We  can  recommend: 

a)  the  analytic  method  -•  formulas  (1)  and  (2)  (for  simplified 
boundary  conditions) : in  case  there  is  no  snow  covering,  there  is  no 


# 


Consequently,  there  is  no  need  to  refine  the  result  obtained 
(we  simply  point  out  that  the  refinement  would  prove  to  be  equal  to 
1 3% ) • 

If,  however,  we  adopt  in  the  same  example:  v = 1.7  m/hr  and  y * 
we  obtain 


*.-*•**  X, 2,0- 0.14 
.Y  = 0,28  - 0.14  = 0.14 


then: 


*'  “uaScw  'l00“  ";0&- 

In  this  case  it  is  necessary  to  utilize  a more  precise  solution. 
First,  according  to  equations  (40)  and  (42)  we  determine: 

ryj 


-i  / • 7 „ o 

P_73  y 2>j  ^ — 15,9  kcal/m  -hr; 


J m-iouo 

't  = 1.73 

•4=» — 1 -22 -f  15,9  0,54  *•  — 13,4  kcal/m-hr. 
According  to  computation  equation  (44)  we  find: 

— 13.4 


vinnn  0.54.15,9-5-13,4  — u.-i  -'4VW  ....  . 

4000 jjTs In  - _ u- - -t.- 0.14=  1/0  hn. 


The  refinement  proved  to  be  equal  to: 


3io^  sc°  = 24<,; 


In  conclusion  it  seems  appropriate  to  us  to  note  that  in  his 
work  the  author  was  provided  by  many  valuable  comments  by  B.  V. 
Proskuryakov  and  N.  N.  Petrunichev. 
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